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We consider the stability problem for a unitary N+1 fermionic model, i.e., a system of N identical 
fermions interacting via zero-range interactions with a different particle, in the case of infinite two- 
body scattering length. We present a slightly more direct and simplihed proof of a recent result 
obtained in [^, where a sufficient stability condition is proved under a suitable assumption on the 
mass ratio. 
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I. INTRODUCTION AND MAIN RESULT 


The study of the quantum mechanical many-body 
problem with pairwise zero-range interactions has re¬ 
ceived a considerable attention in recent years as an effec¬ 
tive model describing the behavior of cold atoms near the 
BEC/BCS crossover (i^a and references therein). The 
correct definition of the model, the occurrence of the Efi¬ 
mov effect and the analysis of the stability problem, i.e., 
the existence of a lower bound for the Hamiltonian, have 
been widely studied both in the physical (i, a 
3, 15 1. [ 10 . Ir^ ) and in the mathematical (see, e.g., 
121 1 lOii literature. It is well known that in the 
two-body case the entire class of Hamiltonians with zero- 
range interaction can be constructed and the spectral 
properties are completely characterized (Q) while, on 
the opposite, for more than two particles an explicit char¬ 
acterization is still lacking. Proceeding by analogy with 
the two-body case, one can construct a physically rea¬ 
sonable class of Hamiltonians usually called Skornyakov 
Ter-Martirosyan (STM) operators. As a matter of fact, 
such operators are symmetric but (in general) are nei¬ 
ther self-adjoint nor bounded from below. This happens, 
for instance, in the case of three identical bosons, where 
it was shown in Q that the STM operator admits self- 
adjoint extensions which can be explicitly constructed 
but they are all unbounded from below and therefore the 
system is unstable. 

Here we are interested in the stability problem in the 
fermionic case, that is when fermions of different species 
interact among themselves. For the most general system 
composed by a mixture of N identical fermions of one 
type (with mass mi) and K identical fermions of a dif¬ 
ferent type (with mass mg), the stability problem for the 
corresponding STM operator is open and some results 
are available only in special cases. For instance, a sys¬ 
tem composed by two identical fermions plus a different 
particle is known to be stable if (and only if) the mass 
ratio 


mi 

a = - 


mo 


( 1 ) 


is smaller than the critical value 13.607 (see, e.g., [1], 0]. 
Further results are available in the 3-1-1 case only (0|). 

In this note we review a result obtained in 0| on the 
stability of a system composed by N, with N > 2, iden¬ 
tical fermions plus a different particle. In particular it is 
shown that stability occurs if 


a < ac{N) 


( 2 ) 


where ac{N) is the solution of the following equation 


A{a,N) :=f(7V-l)(i^f 


v^TT^ 


— arcsin 




= 1 


(3) 

Notice that for each N the function A(q;, N) is increasing, 
goes to infinity for a —>■ 00 and A{0,N) = 0. So there 
is exactly one solution adN) > 0 of ([3]) and moreover 
A{a,N) < 1 for a < adN). We remark that only for 
N = 2 the condition (ED is optimal, i.e., ac(2) = 13.607, 
and therefore the result provides a rigorous proof of what 
is already known in the physical literature. 

For N > 2 the condition is surely not optimal since, as 
it will be clear from the proof, the role of the antisym¬ 
metry is only partially exploited. Nevertheless we believe 
that the result can be of some interest since (ED gives a 
first sufficient stability condition which, apparently, was 
not known before. Some numerical values of adN) are 
listed here 


ad3) = 5.291, 
ac(8) = 1.056, 

ac(9) = 0.823,... (4) 

In particular this means that, in the case of equal masses, 
the system is stable if A < 8. 

Our aim here is to present a slightly more direct and 
simplified proof of the above stability condition, limiting 
ourselves to the simpler but relevant unitary case (see 
below). We hope that such a presentation could help in 
clarifying the points of the proof which must be improved 
to obtain a more satisfactory stability condition. 

Let us consider the formal Hamiltonian for a system 
of N identical fermions with mass mi and a different 
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particle with mass mp (we set h = 1) 

^ ~ ~2^ Si=l —2^^xo+MSi=l <^(^0—Xi) . (5) 

The parameter p. G R is a coupling constant which must 
be properly renormalized in order to give a precise mean¬ 
ing to the expression ([S]). We introduce center of mass 
and relative coordinates 

Xc = (tooXo mi x*^ , yi = xo - Xi (6) 

where M = tuq + Nmi. One has 

if — _3_A rno+mi a 

^ ~ 2M^x„ 2momi Z^i=l 

- wS.<, ^y.-Vyi+MEti<5(y.)- (7) 

Therefore in the center of mass reference frame the sys¬ 
tem is described by the formal Hamiltonian 

H = Ho + flj:t,Siy.) (8) 

where p. — Hq is the free Hamiltonian 


antisymmetric under the exchange of any pair of 
fermions, regular for Xp ^ x^, i = 1,... ,iV. The stan¬ 
dard formulation of the boundary condition satisfied as 
|xo — Xi| —>• 0 is (see, e.g., m) 


(()(xp, . . . ,Xjv) = 


I -1 - - (-1) Q(rpi,x,) 

|xp-Xj| a/ 


0(1) (13) 


where Q is a given function depending on (f> and 

mpXp -I- TOiXj 

^- 

mp -I- mi 

Xj — (xi, . . . , Xj_ 1 , Xj^ 1 ,..., x^v) 


(14) 

(15) 


Notice that in the limiting procedure defining the bound¬ 
ary condition (I13p the vectors rpi and x^ are kept fixed. 
Passing to the relative coordinates Y = (yi,-.-,yw), 
the wave function tp{Y) = ^(xo(Y),...,XAr(Y)) is an 
element of regular for Y ^ H. Moreover, set¬ 

ting Q(Y) = (5(rpi(Y),Xi(Y)), the boundary condition 
m satisfied as |yi| —>■ 0 now reads 


Ho = - Ell Ay. - ^ E.<, Vy. • Vy,. (9) 

with a the mass ratio O- The formal Hamiltonian ® 
can be given a precise meaning as a, possibly self-adjoint, 
operator in the natural Hilbert space Lant(®^^^)i i-e., the 
Hilbert space of antisymmetric, square integrable func¬ 
tions on . More precisely, it is by definition a non 
trivial (self-adjoint) extension of the free Hamiltonian Hq 
restricted to smooth functions vanishing on the set 


^(Y) = 


1 1 

|y*l a 

1 1 

|yi| a 


(-l)*+iQ(Y)+o(l) 
(-l)*+iQ(Y) 


Vi 

\i-i-i ( 


yi=0 

+ 0(1) 


- (- 7) (-i)'’^^C(y*) + ®*(yO • + o(i) (16) 


=Ue{i iv}{Y=(yi,...,yiv)eK3iV|y.=0}. 

ao) 

As we already mentioned, among all possible extensions 
a special role is played by the STM operator, due to the 
fact that it is the natural generalization of the well known 
Hamiltonian with zero-range interaction in the two-body 
case. 

For two (different) particles, extracting the center of 
mass motion and denoting by x the relative coordi¬ 
nate, the domain of the operator consists of functions 
ip G which are regular for x ^ 0 and satisfy the 

following boundary condition as |x| —^ 0 

where q G C depends on ip and a S R has the physical 
meaning of a scattering length. 

Moreover, the Hamiltonian acts as the free Hamiltonian 
for |x| y 0. 

The STM extension Ha in our fermionic N + 1-particle 
system is defined in an analogous way. In the original co¬ 
ordinate system Xp, Xi,..., x^r, and extracting the center 
of mass motion, the domain D{Ha) is made of functions 
(p defined on the set 

{(xp,... ,XAr) G Xc = 0} (12) 


where ^ : R^^“^ —>■ C is an antisymmetric function, 
= yi/|yi| and a G R is the two-body scattering length 
corresponding to the interaction of a fermion with the 
different particle. 

We remark that in (|TO|) (i.e., the boundary condition 
written in the relative coordinates) an extra term ap¬ 
pears, depending on the unit vector namely the di¬ 
rection along which the limit y^ —>■ 0 is taken. As a 
matter of fact, such a term does not contribute to the 
energy of the system (see Section H), yielding the same 
expression one would get if that term was absent (see also 
the Remark at the end of Section H). 

Furthermore, Ha acts as the free Hamiltonian outside the 
set H, i.e., 

{HaiP){Y) = (i?oV')(Y), if Y G R^"^ \ H. (17) 

A special role is played by the parameter-free case of 
infinite scattering length, known as the unitary case. We 
shall denote by H the corresponding STM extension, i.e., 
H :=i7oo. 

The main result discussed in this note is the following 

Theorem. In the unitary case the energy form, i.e., the 
expectation value of the energy, is positive for a < ac{N). 
More precisely, for any ip G D{H), 

{'ip,Hip)>Q, if a < ac{N). 


(18) 
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This in particular implies stability for the unitary + 1 
fermionic model. 

In the next Section we derive a suitable expression for the 
energy form. In Section III we start from such expression 
to explain the steps required to prove our result. In Sec¬ 
tion IV we briefly summarize the content of the paper. 
In Appendix A we collect some technical results useful 
to reformulate the domain and the boundary condition 
characterizing the Hamiltonian. 


The proof of (1^ and (pSjl is postponed in Appendix lAl 
Using the decomposition (HD) and the asymptotic behav¬ 
ior ([22]), the boundary condition (HU) in the unitary case 
can be equivalently written as 

lim w(Y) = (r^)(yj) + ©,(y*) • w*. (27) 

IviKo 

We can now derive the expression for the energy form. 
Taking into account (HZl, the decomposition (USD and 
equation (I^Sl) we have 


II. DERIVATION OF THE ENERGY FORM 

Here we derive a suitable expression for the energy 
form. The key point is to represent the domain D{H) as 
the set of wave functions decomposing as 


= w + Qf 


(19) 


where ru is a smooth function and Qf contains the sin¬ 
gular behavior prescribed in (fTBl) . More precisely, Qf is 
the “potential” produced by the “charge” f distributed 
on the planes {y^ = 0}, i.e.. 


= V; V 


N 


(-ly+i 


TT ^ (2^)t^ 


dK e*'^-^G(K)C(kj), 


( 20 ) 

with K := (ki,..., k^v), / denoting the Fourier trans¬ 
form of / and 


G(K) = 


1 


E* kf + k* • k, ■ 


( 21 ) 


Indeed, it is straightforward to verify that 

(-i)*+'C(y*) 


{Gm) = 


|yi 


-(re)(y.) + o(l), (22) 


for |y^| —>■ 0, where is the inverse Fourier transform 
of 

(h) (k,) = ^^(-i)*+u(k,)^"(ko 


l + Q! 

I 

27r^ 


^ (-1)^+1 |dk,G(K)|(k,) (23) 




and 


L(ki,...,kjv-i) = (24) 

Moreover, it is useful to note that the potential Qf satisfy 
the equation 

{Ho Gf) (Y) = 4^ ^(-l)*+ie(y.) <5(y.) (25) 


in distributional sense and then, in particular, 

{Hogf){Y)=0, iiYGR^^\n. (26) 


{ip, Hip) = liin / dyi .. .| 


£—^0 


'|yil>E 


' |yN-|>e 


dyNip{Y){HoiP){Y) 


= lim 
£—^0 


dyi... j 


|yi|>e 


MywUe 


dy^v {w + gf){Y){How){Y) 


= {w, Hqw) 


-I- lim 
£—^0 


dyi...; 


|yil>E 


' |y;v|>e 


dyN{gf){Y){How){Y). (28) 


In the last integral of we apply Green’s identities. 
Denoting Sf = {y^ G | \yi\ = 1}, for e —>• 0 we have 


/ dyi.../ dyiv(^(Y)(i7o«^)(Y) 

'|yi|>e "'|yAr|>£ 

= V [ dy, [ da{uji) 

IrHN-1) Js? 


AGi) 


2=1 


5|y.l 


|yd=£ 


|y.l=£ 


^(1) 


N 


= dy,C(yi)/ dcr(w,)w|| I -fo(l) 

= dy,^(yi) f dcr(w,) hr^)(y,) 

^ JrHn- 1) Jsf V 

+ 0i(yj) ■ ) + o(l), (29) 


where we have used equation (EH), the asymptotics 
and the boundary condition (071) . Observe that in (0^ 
the term o(l) contains surface integrals going to zero for 
e —>■ 0. Note also that, due to the integration over the 
unit sphere Sf, the contribution of the term 0i(yi) • 
vanishes. Taking into account ((75)) and (1771) . we find 

{ip, Hip) = {w,How) 

+ 4^^(-1)*+7 dk,c(k,)(ro(k.). (30) 

^ jR3(iV-l) 

Inserting in (1301) the explicit expression ((75)) of T^ and 
exploiting the antisymmetry property of the charge f we 
finally obtain 

2N 

{ip,Hip) = {w,How)-\ -$(0, (31) 

TT 

where the quadratic form $ is defined by 

$(C)= ^dK (<l>i(C;K) + (iV-l)$2(C;K)) , (32) 
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with 

K = (ki,..., kAr- 2 ), 

= y'dsL(s,K)|^>,K)|2, (33) 

4>2(C;K)= y'dsdtG(s,t,K)|(s,K)|(t,K), (34) 

Note that in the special case N = 2 the extra variables 
K are absent. 


Remark 

As we have seen in the above computation, the presence 
of the term &i(yi) ■ uj^ in dT51) is irrelevant, in the sense 
that it does not contribute to the energy and therefore 
it can be dropped from the outset. In our opinion, this 
simply means that in the boundary condition the limit 
for |xo — Xi| —>• 0 does not depend on the direction along 
which the limit is taken. In order to stress this fact, m 
should be more correctly written in the following way 

^^dcr(w,)<()(xo,...,XAr) 

— Z —i -+ 0 ( 1 ) (35) 

Xq X^I Q. y 
as |xo — Xi\ —^0. 



A. Change of variables 


We define 


P = s + XI ’ q = t + k,, (36) 

Z H- a ^ ^ rv ^^ 


N-2 


N-2 


2=1 


2 + a 


2=1 


and therefore we obtain 
$ 


1 (C it) = 2n^Jdp + |,(p, k:)P , 

<t.teK)= bpdq-JEilWsA) , (37, 


+ q^ + i^p- q + ^(K) 


where 


,7(p,K)=e(p-T^Ek7,K) 


(38) 


D{K) = 


l+3a 

(l+c<)(l+2c<) 


(Efc? 


2 I 2a 


l+3a ^i<j 


E7<, k,; • k. 


B. Expansion in spherical harmonics 

For any / G L^(IR.^) we consider the expansion 


III. POSITIVITY OF THE ENERGY 

From the above expression (1511) we see that the Theo¬ 
rem is proved if we can show positivity of 4>. Since the 
term $1 is positive, the problem is reduced to show that 

(N - 1 )< 1)2 > -C<^1 

for some constant c < 1. A proof of this fact will be 
given here and for the sake of clarity it will be divided in 
several, but elementary, steps. The strategy will be the 
reduction of the form $2 to one which can be diagonal¬ 
ized. This first requires a suitable change of variables; 
then we exploit the rotational symmetry of $2 to per¬ 
form a partial wave decomposition; once the additional 
degrees of freedom are dropped, the problem reduces to 
bound from below a two-particle energy, which can be 
diagonalized by means of the Fourier transform; to con¬ 
clude the proof it suffices then to go back to the original 
expression and show that, if the condition a < is sat¬ 
isfied, $ is positive. It is worth stressing that at the last 
stage of the proof (see, e.g., (H51) l the fermionic symmetry 
of the charges ^ is totally neglected, in order to diago¬ 
nalize the expression. This is clearly not optimal and an 
improvement of the condition would require a differ¬ 
ent approach. In fact the change of variables itself (see 
(1381) 1. which is the starting point of our analysis, make 
the antisymmetric requirement not apparent and there¬ 
fore should probably be avoided if one wants to track 
down the role of the fermionic antisymmetry. 


00 I 

/(P)=E E (39) 

l—O m——l 

where p = (p, 0p,(j)p) and y)™ denotes the spherical har¬ 
monics of order I, m with Z = 0,1,... and m = —I ,...,/. 
Moreover, we denote by Pi the Legendre polynomial of 
order I explicitly given by 

-'L <«) 

Using the above expansion we derive the following de¬ 
composition of $2 in each subspace of fixed angular mo¬ 
mentum 1: 

QQ ^ poo poo 

^’2(C;K) = 27 r^^ / dp dq rjjZ^{p,'K)r]im{q,'K) 

i=0m=-r° 

, f'dy -_ 

7_1 p2 + q2 z ^pqy + DiK) 

00 I 

=^E E (41) 

/=0 m— — l 

It turns out that (for details see Lemma 3.2]) 

Gi(??z 77 i; K) > 0 , for 1 even, (42) 

0 > Gi{mm;K) > G°{pi^) , for I odd, (43) 
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where G° is defined by 

= 27t j dyP!{y)x 


Jo Jo P + 9 + i+^P^lV 

From (H^ and (H511 we then get 


$2(e;K)> ^ ^ GUmm). 


(45) 


I odd m— — l 


C. Diagonalization 


Let us define 


g^{k) = 


v^. 




(46) 


Then 


Gn9) = 27Tj\yPi{y)x 


X / da;ida ;2 


g3xi^(gXl)g3x2^(g^2) 
g2a:i _|_ g2a;2 _|_ gXi+X2 

g2xi-(gXi)g2x2^(gX2) 


1+a 

ri r ^2xi—f„xi\„2x 

= TT / dyPi(j/) / dxidx2 


/-I 


cosh(a;i - 0 : 2 ) + j^y 


The last integral is a convolution and therefore can be 
diagonalized by means of Fourier transform. Using the 
explicit Fourier transform of the kernel (see, e.g., [8|) we 
find for I odd 


G°(y,„)= dkSi{k)\gj^{k,K)\^. 


(47) 


sinh^/c arcsin 

sin Yf^y) 


(48) 


Si{k) = - . . / dyP;(y)^ 

^ '' sinh(^/cj J cos(, 

and the estimate (03) becomes 

‘i>2(C;K)>^^ |dfc5Kfc)hL(^.K)p. (49) 


I odd m— — l 


D. Bound from below 


We notice that, for any fixed I, Si{k) is an even, G°° 
function of k and limfe_).oo Si{k) = 0. Furthermore for I 
odd we can show that Si{k) is an increasing function of 
I for any fixed k (for details see Lemma 3.5]). Then 

Siik)>S,ik). (50) 

Moreover it is easy to see that 
Si{k) > 5i(0) = 47ri^ arcsin - l] (51) 


where 51(0) < 0. Therefore from we have 

OC I p 

<l>2(e;K)>-|5i(0)|^ ^ dk\gjjk,K)\\ (52) 

1=0 m=-l '' 

which can be rewritten in such a way to reconstruct the 
term $ 1 . Indeed 

/dA|y?„(fc,K)|2= (dxe^^\gi,^{e\-k)\^ 

J JR 

|y,^(p,K)|2 . 

Using this estimate in (15^ we find 
$2(C;K)>-|5 i(0)|^^ X 

00 I pQO - - 

^PP^y-^^^P'^ + D{K)\gim{p,K)\^ 

1=0 m=-l''^ 

= /dPV^lii^p" + ^(K)|y(p,K)p 

= -|^i(0)l2^^^i(^;K) 

= _(7V-l)-iA(a,fV)$i(e;K). 

We are now in position to conclude the proof of the The¬ 
orem. From (I5T]) . (15^ and the inequality above, we get 


2N 

{ip,HiP) > -4>(^) 

TT 


> 


^(l-A(a,iV))ydK$i(^;K), (53) 
and taking a < adN) we obtain the desired result 

{d,Hd) > 0. 


IV. CONCLUSIONS 

We have reported on a derivation of a sufficient con¬ 
dition on the mass ratio for the stability for the unitary 
fV -|- 1 fermionic model. Such a condition, which is op¬ 
timal only in the two-particle case, is nevertheless non 
trivial for generic N. For istance it provides stability in 
the case of equal masses up to A^ = 8. We have also 
described the main steps of the proof, enlightening the 
points to be improved to get a more refined stability con¬ 
dition. 


Appendix A: Properties of the potential 


We first prove the asymptotic expression (12^ for the 
potential for jy^j —>■ 0. It is convenient to isolate the 
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z-th term of the sum in (1^0)1 


A : = 


1^2^ Jdk,e^y-^^G{K) 


dk 


g*yi-k. 


f e»- -.(k. a + l,) \ 

y '111(111 +ki a+ li)j’* * 


where 


a= 

^ = k? + T,n<l,n k, • ki. 


T+K l^j<l-,jMi 
By an explicit computation we find 


/dk 


giy.-k. 27r2 


and as |yi| —>■ 0 (see [HI) 

giyi-ki^jj. . _|_ ^^ 

kf (kf + kj • a + b) 


/dk 


k? 


= /dk 


(A2) 

(A3) 

(A4) 


kj • a + 6 


kf (kf + ki • a + b) 

/•„, (e-y-->---l)(k,-a + 5) 

J ' k2(k2 + ki-a + 6) 


= TT^ ^746 — a^ + o(l) 
Using (EH) and (ESI) in EB) 




Considering the contribution of the other terms of the 
sum in (ESI), it is now easy to derive the asymptotic be¬ 
havior ([22]). 

Let us verify that the potential satisfies equation 
(ESj)- From the definition (1201) and the expression of Hq 
in Fourier transform we have 




(2^) 


%N 


N 


E(- 

i=i 


R3W 


dKe'^'^^(k 


N 


/ dk,"*' 


■Vj 


(2^)2 V 

= 47r^(-l)*+iC(y0%0 


(A6) 


A = 


2(-l) 


1+1 


TT (27r)2^ 
27r2 


dk, e*y-‘^-7(k,)x 

L(ki) + o(l) 


2 ^/1 A 2a 


1 1 ~ 27r' 

|y»l 1 + a 
(-i)*+'7(y.) VT+^ (- 1)*+1 


|yi| l + a (27r)5(^-i) 

X /dk,e"y--^-7(k,)L(k,) + o(l). (A5) 
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